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Abstract 



In an earlier work, Damle and the author (Phys. Rev. B in press; cond 
mat/9705206| ) demonstrated the central role played by incoherent, inelastic 
processes in transport near two-dimensional quantum critical points. This 
paper extends these results to the case of a quantum transition in an anyon 
gas between a fractional quantized Hall state and an insulator, induced by 
varying the strength of an external periodic potential. We use the quantum 
field theory for this transition introduced by Chen, Fisher and Wu (Phys. Rev. 
B 48, 13749 (1993)). The longitudinal and Hall conductivities at the critical 
point are both e'^/h times non-trivial, fully universal functions of huj/kBT (w is 
the measuring frequency) . These functions are computed using a combination 
of perturbation theory on the Kubo formula, and the solution of a quantum 
Boltzmann equation for the anyonic quasiparticles and quasiholes. The results 
include the values of the d.c. conductivities (fao/kBT 0); earlier work had 
been restricted strictly to T = 0, and had therefore computed only the high 
frequency a.c. conductivities with hu/kBT oo. 
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Typeset using REVT^ 



I. INTRODUCTION 



The remarkably rich zero temperature phase diagram of a two-dimensional electron gas 
in a high magnetic field offers an ideal laboratory for exploring the physics of systems 
near a quantum critical point. Much theoretical and experimental attention has focussed 
on the transitions between the various quantum Hall plateau states and the insulator 0J^. 

The vast majority of the theoretical work has studied transitions between integer quan- 
tum Hall plateaus using a model of non-interacting electrons undergoing phase-coherent 
transport in the presence of a random potential 0,^: this transport can be described by 
non-linear sigma model field theories |^ or by a network model 0. In this scenario, interac- 
tions only provide a phase-breaking length which cuts off the critical quantum interference 
effects, but are irrelevant perturbations at the quantum-critical point. It has also been 
suggested that the transitions involving fractional quantum Hall plateaus can be described 
using closely related models of non- interacting composite fermions [Tj-Q. 

However, there have been a number of indications in recent experiments that interactions 
play a more fundamental role 0|l^JTl[]. Measurements of the non-linear transport ll2|,^ have 



allowed direction measurement of the dynamic critical exponent z\ the measured value 2 ~ 1 
is incompatible with critical points of non-interacting electrons, but can possibly be under- 
stood to be a consequence of the Coulomb repulsion at an interacting critical point |T^|T^. A 
more direct indication of the importance comes from measurements of frequency dependent 
transport. It was found that the width of transition region between the plateaus scales with 
the same power of frequency, cj, and temperature, T |T^: such a property only holds at quan- 
tum transitions described by a fixed point with non- zero interactions, while non-interacting 
fixed points with irrelevant interactions have very different sensitivities to perturbations in 



ijj and T ||T5|. 



Non-zero temperature transport near interacting quantum critical points in two dimen- 
sions was studied by Damle and author in a previous paper |Tl| (hereafter referred to as I). 
A number of general observations were made, supported by an explicit calculation on the 
superfluid-insulator transition in a Mott-Hubbard model for bosons. In particular, it was 
argued that the conductivity at the critical point was given by jh times a non-trivial uni- 
versal function of the ratio, huo/ksT. Attention was drawn to the rather different physical 
properties of the limiting regimes of frequency and temperature: (i) fiu) ^ ksT with phase 
coherent transport of excitations created by the external electric field, and [ii) Tiuj <C ksT, 
with incoherent, quantum relaxational transport involving repeated inelastic scattering be- 
tween pre-existing thermally, excited carriers, and relaxation to local equilibrium. The d.c. 
transport is clearly in the latter regime, and it was noted in I that the d.c. conductivity was 
universal even though it was dominated by inelastic processes: this was possible because the 
cross-section for scattering between the carriers became universal at the quantum critical 
point. 

This paper will extend the considerations of I to one of the simplest possible interacting 
quantum critical points in a fractional quantum Hall system. Any experimentally realistic 
model must account for both the Coulomb interactions among the electrons and for the ex- 
ternal random potential. Such a model is of forbidding complexity, and we shall choose here 
to consider non-random models in the presence of an external periodic potential with only 
a statistical gauge force interaction between the excitations. Our choice is motivated by the 
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desire to explore seriously the consequences of interactions in the simplest possible context, 
and the arguments above that inelastic, incoherent processes are essential to any analysis of 
the experimentally observed critical properties. The absence of a random potential means 
that there are no localized states: these are essential for the appearance of quantum Hall 
plateaus as a function of the external magnetic field, as they appear when the Fermi level 
is in a region of localized states. However, the localized states are not expected to be as 
crucial in the experimentally relevant critical region (the studies of non-interacting electrons 
in a random potential examine the complicated physics of the divergence of the localization 
length at the critical point in perfectly phase coherent transport, but this is always cut off in 
interacting models at some finite scale due to incoherent, inelastic processes present at any 
non-zero temperature). A second motivation for examining the case of the periodic potential 
comes from its importance in the early studies of the integer quantum Hall effect |T6[: the 



topological nature of the quantization of the Hall conductivity first emerged in this context. 
A number of earlier studies have examined the fractional quantum Hall effect in a periodic 



potential. Kol and Read [|T^ examined general global criteria that must be satisfied by such 
states. Wen and Wu |]18[ and Chen, Fisher and Wu first examined the quantum critical 
points between different fractional quantum Hall states and the insulator, and developed a 
low energy field-theoretic descriptions of the transitions. They argued that these transitions 
could be understood in terms of a Mott transition of an anyon gas in a periodic potential 
and an external magnetic field. The field theory of the anyon gas consisted of a U{1) 
Chern-Simons gauge field coupled to either relativistic bosonic or Dirac Fermi fields. In 



an interesting recent work, Pfannkuche and MacDonald pO| have numerically shown the 
presence of a transition between a fractional quantum Hall state and an insulator induced 
by varying the strength of an external periodic potential acting on an electron gas in a 
strong magnetic field: they also present evidence for scaling behavior in the spectrum of 
small systems near the critical point. 

This paper shall describe the nature of non-zero T and u transport in the U{1) Chern- 
Simons gauge field coupled to Dirac Fermi field model of Chen, Fisher and Wu [|1^]. Our 
motivation for chosing this particular model is mainly technical: it provides the simplest 
possible context for describing the physics of incoherent transport in quantum Hall critical 
points, and allows us to keep technical details to a minimum. Our results can also be 
extended to models of bosons coupled to a Chern-Simons gauge field using a expansion, 
and this will be discussed in Appendix 0. We will show that at the critical point, the 
longitudinal (a^x) and Hall {(Jxy) conductivities obey the scaling forms 



= ; = (1.1) 





We shall extend the quantum transport theory of I to obtain explicit expressions for the 
universal functions Tj^x and E^.^. The earlier computations of these conductivities were 
carried out precisely at T = 0, and therefore yielded the values of ^^^^.(cxd) and J2xy{oo) which 
describe phase-coherent transport. The d.c. conductivities are however given by S^a,(0) and 
Sa;j/(0): these are quite different, but also universal numbers, characterizing incoherent, 
relaxational transport, and will be computed in this paper. 

As many of the considerations of this paper are similar to those of I, we now explicitly 
point out some aspects which are different. The central focus of our analysis shall not be 
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the conductivities a, but related irreducible conductivities a\ these are diagrammatically 
defined as being irreducible towards cutting any single Chern-Simons gauge field propagator. 
We shall find that many of the issues arising in the computation of similar to those 

discussed in I: the focus shall be on a quasiparticle-quasihole contribution which is dominated 
by collisions among the excitations, and is best described by a quantum transport equation. 
The physics going into the computation of axy is quite different. The real part of axy is non- 
dissipative and is actually best computed by a simple perturbation theory, which already 
yields a result which is an interesting function of huj/kBT. Indeed, such a perturbation 
theory was presented in Ref at T = 0, and here we consider the orders of limit of a; — >■ 
and T ^ more carefully. Our final results for a^x also has a form rather different from 
that obtained in I: it is given by complicated interplay between different contributions to 
a: in particular we find that in the low frequency collision-dominated regime it is possible 
to have a^x be an increasing, rather than decreasing, function of hu/kBT. For the latter 
situation we find that axx{^) has a peak at a; ~ ksT/h, and eventually decreases for large 
enough u. 

As noted above, an essential ingredient in our computation is the quantum transport 
equation for quasiparticles and quasiholes. These excitations are anyons, and so their scat- 
tering cross-section is singular in the forward direction: for momentum transfer q, the cross- 
section ~ |2T|. This scattering cross-section appears in the collision term of the trans- 
port equation, and could potentially cause infrared divergences. However we shall show 
explicitly that the potential divergences all cancel, and the physical transport co-efficients 
of the anyon gas are well defined. This is an important, albeit technical, consistency check 
of our analysis, and occupies a significant part of Section fV[ 

The outline of the remainder of the paper is as follows. We will introduce the model and 
discuss some important general issues in Section y. The perturbation theory for a will be 
presented in Section |I|: this will allow us to identify the different physical contributions 
to a. One of these, the quasiparticle/quasihole contribution is a delta function at zero 
frequency at leading order in a. This is an artifact, and has to be repaired by an infinite 
order resummation: the latter is most conveniently carried out using a quantum transport 
equation analysis which will be discussed in Section 0. We will synthesize our results for 
a into those for the physically observable a in Section 0, and conclude with some general 
remarks in Section |V|. 



II. THE MODEL AND SOME GENERAL CONSIDERATIONS 



To the extent possible, we will follow the notation of Ref |r^|. We will be working the 
relativistic models with a velocity c, and use units with h = c = ks = ^, except where we 
explicitly quote conductivities in units of e'^/h. The imaginary time action of Dirac fermions 
coupled to a f/(l) Chern-Simons gauge field introduced in Ref |]19 



IS 



l/T 



dr I dx<j tp-^ [7^ {d^ - iga^ - iqeA^) + iMi] tpi 

1 



(2.1) 
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where x = {x, y), and the Greek indices extend over x, y, t. The Chern-Simons gauge field is 
a^, the external vector potential driving the currents in A^, and 'i/'i^2 are two species of Dirac 
fermions with charge qe. The Dirac matrices are 7^ = i{—Ty,Tx,Tz) where r are the Pauli 
matrices, and ip = V'Wt- Notice that we are using the symbol r to indicate both imaginary 
time and the Pauli matrices, but the meaning should be clear from the context. The sole 
role of and the Chern-Simons term is to attach flux tubes onto each Dirac fermion and 
to convert it to a Dirac anyon with statistics parameter (1 — a) where 

we are using the nomenclature where a statistics parameter of 1 corresponds to fermions, 
and to bosons. Notice that we have two species of anyons, ipi^2'- these are required by 
the doubling theorems for chiral fermions on a lattice. The anyons/anti- anyons represent 
the Laughlin quasiparticles/quasiholes of the fractional quantized Hall state. So to describe 
physics in the vicinity of the simplest Laughlin quantum Hall state at flUing z/ = 1/3 we 
should choose a = 2/3 and q = 1/3. As we will review below, the action S describes |T9| 
a quantum transition between fractional quantized Hall state with T = conductivity (the 
Latin indices extend over x, y) 



and an insulator with aij = 0. Notice that, as expected, (|2.3|) gives a Hall conductance of 
eV(3/i) for q = 1/3 and a = 2/3. 

We now make some general remarks about the nature of charge transport in Chern- 



Simons theories: these results are well known [p|,p^ , but our formulation is especially suited 
to our subsequent analysis of inelastic processes using the quantum Boltzmann equation. 

From the equation of motion of implied by S, we see that role of is to impose the 
constraint 

(^t^ux^uax = g (V'i7m^i + V'27mV'2) • (2.4) 

This constraint will be imposed exactly in all stages of our analysis. The r component of 
this constraint is 

eijOittj = -g (V^iV^i + ^2^2) • (2.5) 
where the Latin indices extend over x, y. The spatial components of ( |2.4D can be written as 

- qetijEaj = g {Jii + J2i) (2.6) 
where Eq is the 'electric fleld' associated with the gauge field 

Eai = drai - dittr, (2.7) 

and Ji, J2 are the electrical currents carried by the two species of Dirac fermions: 
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(2.8) 



and similarly for J2. 

We shall be interested in the response of S to an external physical electric field E given, 
of course, by 

= d^Ai - diA^, (2.9) 

The currents induced by such a field, in turn, induce a non-zero value of the E^ field, by 
(p.6|). It is therefore useful to define 

Et = E+— E, (2.10) 

qe 

as the total effective field experienced by the Dirac fermions. The field Ey will appear 
naturally in our subsequent formulation of the quantum transport equations. It is then 
useful to consider the response of the fermions by introducing the linear response relations 

Jii{uj) = aiij{uj)ETj{uj) (2.11) 

and similarly for J 2- We have also introduced a frequency u to allow for the response to a 
time-dependent electric field. In diagrammatic terms, if we compute a^. via a Kubo formula, 
we should only include diagrams which are irreducible towards cutting any propagator: 
the irreducibility requirement is equivalent to considering the response to the effective field 
Et^. 

The experimentally measured response, aij, of course, relates the total current to the 
external electric field 



Jii{uj) + J2i{uj) = aij{uj)Ej{uj) 



(2.12) 



From ( |2 . 6| , |2 . 1 0| , |2 . 1 l| , |2n^) we can obtain the following basic relationship P,[T9| between the 
response functions 



-In 



(2.13) 



When we are well away from the critical point, |Mi_2| is large, and the ai^2 are easily 



computed at T = 0; their exact values are [|I9 



eijSgn(M 



1,2; 



(2.14) 



We see therefore that with the choices Mi > and M2 > 0, ( p.l3| , prT4D lead to the con- 
ductivity (|2.3| ) of the fractional quantized Hall state. The system undergoes a transition 
to an insulator if Mi changes sign while M2 remains positive: then the total conductivity 
vanishes. 

The remainder of the paper will study the critical point where M2 is large and positive 
(and therefore the ip2 fermionic excitations have large gap) while the renormalized Mi = 
at T = 0. The 5^2 conductivity therefore continues to satisfy ( |2.14| ) at all low T and a;; all 
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corrections are suppressed by factors of e"^^/^. The response of the ipi fermions is much 
more non-trivial, and will occupy the remainder of the paper to elucidate: it has a highly 
non-trivial, universal function of u/T. As our attention will focus solely on the response of 
the ipi fermions, we will henceforth drop the subscript 1, and simply refer to them as the ijj 
fermions. 

The central purpose of this paper is to analyze the quantum transport properties of the 



model S in ( |2.1| ) at T > 0. While the action S is Lorentz invariant, the presence of a 
nonzero T breaks the Lorentz symmetry, and the space and time dependence of observable 
correlators will be quite different. Therefore, there is no particular advantage to working 
with a covariant gauge, and it pays to choose a gauge better suited to physics at T > 0. It 
turns out that by far the most convenient choice is the Coulomb gauge: the perturbative 
calculations to be considered in Section |ITT| are the most straightforward, and the quantum 
transport equations of Section |rV| are readily derived in this gauge. The Coulomb gauge 
corresponds to the choice 

d,a, = 0, (2.15) 



combined with an exact integral of S over a^. Then the constraint ( |2.5| ) can be solved to 
give 

a. = <7^^V (2.16) 

where d"^ = didi. Note that this is an equal-time relationship-the absence of retardation in 
the gauge field interactions is an important simplifying feature of this gauge. Inserting ( |2.16| ) 
into (|2.1| ) (and dropping the non-critical il)2 fermions) we obtain the Hamiltonian density 



H = Ho + Hi 



Hi = j 



\ljH-irA-M^T,)i?\ (2.17) 



(2.18) 



We have inserted a bare mass Mq in to ensure that the system is at its critical point, 
the value of Mq will have to be adjusted order by order in a so that the renormalized mass 
vanishes at T = 0. Note that the gauge field has disappeared as an independent dynamical 
degree of freedom: each Dirac particle/anti-particle now has a flux tube attached, and it is 
this flux tube which is responsible for the long-range force in Hi. This economical description 
of the degrees of freedom will be quite useful when we consider the transport equation at 
nonzero temperature. 

The remainder of the computations in this paper will be carried out using the Hamilto- 
nian H. 



III. PERTURBATION THEORY AND KUBO FORMULA 



We will begin in Section [III A| , as in the analysis for the bosonic case in I and in Ref p2 



by an examination of the perturbative structure of the single particle self energy. Then, 
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we will proceed in Sections [IIIB| and [III Q to an evaluation of axx and Jxy using the Kubo 
formula. We will gain an understanding of the physically different contributions to a by 
examining the structure of the bare perturbation theory in powers of a. The results for a 
will have some unphysical features which will finally be rectified by an analysis using the 
quantum Boltzmann equation in the following Section P^ . 



A. Self energy 

The purpose of this section is to understand the structure of the mass renormalization 
in a single fermion propagator. As noted earlier, we need to adjust the value of Mq so 
that the renormalized mass vanishes at T = 0. However the renormalized mass need not 
vanish at T > 0, and we will describe its universal T dependence here. We will denote the 
renormalized mass by M(T), and clearly we should have M(0) = 0. 

To leading order in a, it turns out to be sufficient to consider evaluation of M{T) in a 
bare perturbation theory in a, without the need for any self-consistent mass renormalization; 
the latter only modify results at higher orders in a. As Mq is expected to be of order a, we 
may therefore set Mq = within the propagators that appear in the self energy. To first 
order in a, the fermion self energy for H is 

S(k) = 27rme,,T V / -^^^^Tr ( -r, - r, -] . (3.1) 

The trace is over the Dirac spinor space, e„ is a fermionic Matsubara frequency, and to this 
order, E is independent of the fermion frequency. Evaluating the trace and performing the 
frequency summation, we get 

E(k)^-™r./|^^^tanhO,/2T) (3.2) 

The leading shows that S(0) is precisely a mass renormalization, and therefore 

M(T) = Mo+S(0)|^^o. (3.3) 

Using M(0) = 0, we get 

where A is an ultraviolet momentum cutoff, whose nature we will not have to precisely 
specify. Finally, inserting (|3^ ) and (|3.2|) into ( p.3|) we get 



where f^{p) is the Fermi function 
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The integral in ( p.5|) is easily evaluated, and we get our final result 

M(T) = -2aTln(2). (3.7) 

We will find later that in some cases it is necessary to include this renormalized mass in 
the fermion propagators, while in other cases we will simply be able to use the zeroth order 
result M{T) = 0. 



The leading contribution to a^x is given by a single fermion polarization bubble, with the 
mass M{T) above included in the propagators. To first order in a there are three additional 
diagrams, considered in Ref I1T3, which contribute to the conductivity; however, these only 
yield contributions to axy and will be discussed in Section |1II Q . The fermion polarization 
contribution to arr is 



2^2 



[lUOr. 



W) 



rTr r^- 



-2e„ + k-f-M(T)r,) 



1 



(3.8) 



;-2(e„ + c^0 + k-r-M(T)r,); 
The trace is over the Dirac spinor space, and evaluating it and simplifying a bit, we obtain 



2„2 



(Txx\li^r: 



27cq^e 

hUr, 



where 



Akl + ul 



(3.9) 



el = k^ + M\T) (3.10) 

The two terms in in the integrand in ( |3.9| ) cancel at T = in the ci;„ — *• limit, but this is 
not evident in the present form. To make this explicit, insert 1 = dkx/dkx in front of the 
second term, and integrate by parts under the kx integral; this yields 



2^2 



hUn 



dpk 



(2vr)M4 + 4; 



Akl + ujI 



^1 + 



(3.11) 



We can now evaluate the summation over and analytically continue the result to real uj; 
as in I it is convenient to separate the result into two distinct pieces 



0", 



a 



(3.12) 



The first piece, a'^P{u), is the contribution of thermally excited quasiparticles: 



a. 



qp 



(a;) = S., 



2„2 r 



h 



tt6{uj) 



iV { - 



(Pk ki ( dfiskY 



(2vr: 



,2 c2 



dei 



Ti5{uo/T) + iV 



T 



(3.13) 
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where A^(a), obtained after using (|3.7|), is a pure number 



POO 

Af{a) = dy 

J 2a In 2 



Ay cosh^(y/2) 

= \n{2) + 0{a^) (3.14) 

At this order, this contribution to the real part of the conductivity is simply a delta function 
of frequency of weight oc T (which satisfies the scaling form (|1 . 1| ) ) . A more sophisticated 
quantum transport analysis (of the type carried out in I), to be discussed in the next section, 
will show how this delta function is broadened out by inelastic, incoherent collisions that 
appear at higher orders in a. The second piece, a'^°^{uj) is a smooth continuum contribution 
from the creation of quasi-particle/quasi-hole pairs by the external source, to the real part 
of the conductivity, and contains the coherent contribution in the limit u/T —>■ oo; after a 
change of variables in the momentum integration from k to Ek, it can be written as 

~coh. . Q'e' foo -zu;{el + M\T))tanh{e,/2T) 



2h J\M{T)\ ^fc(4£:| - {uj + ir]f 

where ?7 is a positive infinitesimal. The real part of this result can be evaluated in closed 
form 



Re 



'^:fi^)] = ^^(1^1 - 2M(T)) [l + tanh(|a;|/4T), (3.16) 



while the imaginary part can be obtained by a Kramers-Kronig transform of Re 
7rg^e^/8/i. Note that, after using ( |3.7| ) for M(T), the result ( p.l5| ) is consistent with the 



scaling form (|1.1| ). We show a plot of the frequency dependence of the real and imaginary 
parts of a™^ in Fig|l|. There are weak threshold singularities at a; = 2M(T): a discontinuity 
in Re and a corresponding logarithmic divergence in Im[aa;a;]. These are artifacts of 
the absence of damping in the fermion propagators at this order, and the singularities are 
expected to be smoothed out at higher orders. 

The previous computation of Re [axx] in Ref was carried out at T = 0, in which we 



obtain only the uj/T ^ oo limit of cr™^ which, from ( |3.15| ), is given by 



a^^{iu/T ^ oo) = nq'eySh. (3.17) 

In the limit u/T relevant to d.c. transport, a^x is dominated by a^^, which at this order 
is a delta function at zero frequency: as stated earlier, this will be repaired in Section |I^ by 
a quantum transport analysis. 



C. (7xy 

Unlike a^x, the zeroth order in a result for a^y vanishes identically at the critical point. It 
is necessary to go to first order in a, where one finds two distinct types of contributions. The 
first (denoted d'^jj) comes from the simple fermion polarization bubble considered above for 
axx, but with the first order result for the mass M{T) in ( |3.7|) included. The second comes 
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from the diagrams, considered in Ref. [0, which are exphcitly first order in a: these account 
for the momentum dependence of the self-energy (denoted cr^^^) and the corresponding vertex 
correction (denoted o. 



(3)> 
xy 



Considering first the contribution of the simple polarization bubble. This gives 



2„2 



huJr, 



■Ti: 



rTr r^- 



1 



-ze„ + k-f-M(T)r, 



X r„ 



(-^(e„ + + k ■ f - M(T)r,) ) ' ^^'^^^ 

Evaluating the Dirac trace and the frequency summation, we find that the leading 1/uJn 
explicitly cancels out: as a result there will be no singular 5{uj) contributions to a^y as 
there were to a^x- We change variables in the momentum integration from k to e^, and 
analytically continue to real frequencies, and the answer takes the final form 

2q\^M{T) foo tanh(efe/2T) 

h J\MiT)\Uel-{u + tvy ^ ' 

Notice that there is M{T) in the prefactor, and so this result is at least first order in a. 
Nevertheless, it is not permissible to simply set M[T) = in the lower limit of the integra- 
tion to obtain the leading result: this is because the integrand has a potential logarithmic 
divergence for ujn = 0, which is cutoff by the lower limit. The result (|3.19| ) vanishes as 
~ iT /uj in the limit uj/T oo, while in the d.c. limit uj /T ^ we get after using ( |3.7|) and 
manipulations described in Appendix ^ 



a 



{uj/T ^ 0) 



a ln(l/a) 



ln(2) , 



37-7 ln(2) - 3 ln(ln(2)) - 36C'(-1)} a + 0{a^ 



(3.20) 



Here 7 = 0.5772... is Euler's constant and C,{s) is the Reimann zeta function. Notice 
the non-analytic aln(l/a) term: the proper T dependence of the mass M{T) discussed in 
Section [111 A| was required to obtain this. 

We now consider the three diagrams for a which are explicitly first order in a. These 
were considered in Ref [|19[ at T = 0; the T > calculations are forbidingly complicated in 



the Lorentz gauge used in Ref [0 because of the retardation in the gauge field propagator; 
they simplify considerably in the Coulomb gauge we are using here. To first order in a 
we can set M{T) = in the propagators of these diagrams; this is justified a posteriori 
by the absence of infrared divergences in our evaluation of these diagrams. The first such 
contribution is the sum of the two diagrams accounting for the momentum dependence of 
the self-energy correction in one of the fermion propagators: 

^xyi^^n) = 




(3.21) 
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The second term in the first square bracket is the subtraction that corrects for the 
momentum-independent mass renormalization that was aheady accounted for by including a 
M(T) in . The final contribution to a^y is the first-order vertex correction corresponding 
to the self-energy in M"^^- 



f(3) 



xy ■ 



An^q^e^a . 2 /" " ^ " P i^i ~ Pi) 



X Tr 



■ 1 [ 1 1 

{-i€n + k ■ r) \ {-iQn + p-r) " (-^(^^„ + + p ■ f) ^ 

1 ill 



iil„. + p ■ r) {-i{iln + ujn) + p ■ t) j (-i(e„ + k ■ f) 



(3.22) 



We now evaluate the trace over Dirac indices, sum over frequencies in ( p.21|) and analytically 
continue to real frequencies to obtain: 

^(2) 87rVe^" f d^k (Pp tanh(p/2T) tanh(fc/2T) /^(p^ _ k • p) 



' h J (27r)2(27r)2 pk{Ak^ - {uj + ir^Y) \ (k - p)^ ) 

2,.e^„ r r ^M2T)i.mki^)^ ^^^^^^ 



h Jo Jp (4/^2 — {uj + irjy 

where in the second expression we have evaluated the angular integrals. As discussed in 
Ref |]TU[, it is important in ( |3.23D that the k integral be evaluated first as it is convergent. 
There is a logarithmic divergence in the subsequent p integral but this will cancel against 
an opposite contribution in aj^J. The trace and summation over frequencies in (|3.22|) yields 

,(3) ^ 327rVe^" f (Pk (Pp (k x p)^ tanh(p/2T) tanh(A;/2T) 



' h J (27r)2(27r)2pA;(k-p)2(4P-(cj + ir7)2)(V-(cj + ir7)2) 

_ _4gV« .00 .00 tanh(p/2T) tanh(fc/2T) 

- h Jo "^^Jo ^^^^^^^^^ '^\4P-(a; + ^r/)2)(4p2-(^ + ,,^)2) ^^'^^^ 

where again the second expression is obtained after evaluating the angular integrals, and we 
choose to evaluate the k integral first. The integrals can be evaluated explicitly in the limit 

u/T 00, and we get 



In the opposite limit, u/T ^ 0, we evaluated the integrals numerically, and found to a very 
high accuracy that 

- 0) + ag)(../T ^ 0) = (3.26) 

although we have no analytic proof of this simple result. The full frequency dependence is 
most easily evaluated by first taking the imaginary parts of (|3.21|) and ( |3.22|) , which can be 
shown to simplify to 
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Im 



— a;tanh(|cu|/4T) / dpf{p/T)V 
h Jo 



4p2 — cu^ 



(3.27) 



The real part is obtained by adding the constant g^e^a(7r^ + 4)/16/i to the Kramers-Kronig 
transform of ( p. 271 ). 

We now combine ( |3.19|) , ( |3.23|) and (|3.24|) to obtain the final results for a^y We have 
the limiting values 



axy{u/T 0) 
a^y{u/T oo) 





[(■f) 


h 






Y7r2 + 4 


h 


A 16 



a ln(l/a) + (0.650988 . . .)a + 0{a^) 
) a + 0(a=^) 



(3.28) 



The second of these results has been obtained previously [lT9| , p3| ; that calculation was done 
in a Lorentz invariant gauge, and it is reassuring that the same result has been obtained 
here in a completely different, non-relativistic gauge. Notice that there are no delta function 
contributions to (Jxy'- this implies that the entire contribution is due to coherent quasi- 
particle/quasi-hole creation, and the coUisionless flow of pre-existing thermally excited par- 
ticles/holes contributes only to a^P 

We conclude this section be showing a plot of the frequency dependence of the real and 
imaginary parts of a^y in Fig ^ Notice that, as in axx there are threshold singularities 
at a; = 2M(T): a discontinuity in Im [aa-a,] and a corresponding logarithmic divergence in 
Im \axx\- Again these artifacts of the absence of damping in the fermion propagators at this 
order. 



IV. QUANTUM TRANSPORT EQUATIONS 

The perturbative analysis of a in Section |ITT| produced results which had singularities at 
two frequencies: a delta function in at u; = 0, and much weaker logarithmic singularities 
at w = 2M(T). Both these features are in fact artifacts of the low order results, and 
disappear when damping of fermion excitations is accounted for at higher orders. In this 
section, we shall show how the delta function singularity at = is broadened out at order 
a^. The singularity at u; = 2M(T) is due to a threshold to creation of particle-hole pairs: 
it is also expected to be smoothed out at order a^, but determination of the precise form of 
this is left as an open problem. 

So to reiterate, this section will only describe the changes to the result for a^P in ( ^.13[ ) 



at higher order in a. The perturbative results for the other pieces of cr, in particular, (|3.15|) 



for axx and (|3.19| ), (|3.23|) , and ( p. 24]) for a^y, will be used unchanged in our subsequent 



considerations as their spurious singularities are relatively minor. 

The delta function singularity in cf^P, will be explicitly shown below to be due to the 
coUisionless transport of thermally excited quasiparticles and quasiholes. This will be done as 
in I: by writing down transport equations for the distribution functions of these excitations. 
The coUisionless form of these equations appear in Section p.V A| and their solution reproduces 



the result for ct^p obtained in Section [IIIB| ; then, in Section |VB|, we include collision terms 
and present the numerical solution. 



13 



A. Collisionless transport of quasiparticles and quasiholes 



We begin by recalling the defining relation for a: it is the linear response in the 

current to the total internal field E^. Fortunately, the field appears in the transport 
equation in a simple and natural way: the effect of the Hartree and Fock terms is precisely 



such as to modify the external field E to Ey As a result we can simply omit these 

terms while replacing E by E7-, and determination of the a becomes, as in I, a matter of 
balancing the ballistic, particle streaming terms with the collision term. 

Throughout Section |^ we will work with fermionic excitations with a spectrum Sk = k 
i.e. we will set the mass M{T) = 0; as the Hartree and Fock terms in the self-energy are also 
implicitly accounted for in our analysis, this is equivalent to simply disregarding such terms 
and working with a Hq with Mq = 0. In this approximation the transport analysis will 
yield a normalization constant for the spectral weight in the conductivity in ( p.l3| ) which 
is given by the result correct to leading order in a, M{a) = ln(2). Our objective will then 
be to describe the broadening of the delta function in ( |3.13| ), while omiting all terms which 
lead to corrections in its total spectral weight. 

The simplest formulation of the transport equations is in a basis which diagonalizes the 
Hamiltonian Hq. To do this, we first express ip in its Fourier components 

i'(^A-f^J':'^A)e-^', (4.1) 



(27r)2 Vc2(k,t) 

and then perform a unitary transformation from the Fourier mode operators (ci,C2) to 
(7+, 7-): 

Ci(A:) = ^(74k)+7-(k)) 

C2(A:) = ^(7+(k)-7-(k)). (4.2) 

We have introduced here a notational convention that we shall find quite useful in the 
following: as k is a two-dimensional momentum, we can define the complex number K by 

K = kj. + iky where k = {k^., ky) (4-3) 

and = |k| = \K\. Expressing the Hamiltonian Hq in terms of 7-1-, we obtain the simple 
result 

^0 = E / (0^^7l(k)7A(k) (4.4) 

where the sum over A extends over +, — . So 7+ creates a particle with energy A;, while 7^ 
creates a particle with energy —k. At zero temperature the negative energy states will be 
filled: so we identify 7+ as the creation operator for quasiparticles, while 7_ is the creation 
operator for quasiholes. Notice that the quasiparticles and quasiholes are distinct excitations 
in different bands, and the absence of a quasiparticle is not equivalent to the presence of a 
quasihole: this latter feature was also present in the bosonic analysis of I, and as argued 
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there, allows the system to acquire a non-zero d.c. resistance while maintaining conservation 
of total momentum. 

Let us also, for future use, express the interaction Hamiltonian Hi in terms of the 7±: 



^1= E 

A1A2A3A4 



cPki (fk2 (fkg 
(27r)2 (27r)2 (27r)2 

X TA,A2A3A4(ki,k2,q)7l4(ki + q)7l3(k2 - q)7A2(k2)7Ai(ki) (4.5) 



where 



7AiA2A3A4(ki, k2, q) — 



i + A1A4— - — ■ — - — 
|ki + q\ki 



Q{K; - Q*) Q*K2 
A3 — n i A2 — ; — 



(4.6) 



Notice that the scattering amplitude, T, is singular for forward scattering: T ~ This 
is a characteristic property of scattering between anyons |^T|]. We will have to keep careful 
track of this potential infrared singularity in our subsequent transport analysis. 

Finally, we also express the electrical current, defined in ( |2.8| ), in terms of the 7±. For the 
case of a spatially independent current (which is the only case of interest here), the result 
can be written as 

J = J/ + Jii (4.7) 

with 

d^k Ak 



(27r)2 k 



■7i(k)7A(k) (4.8) 



and 

d^k (z X k) 



fii = -iqe 



J±±_y±±^ (7l(k)7-(k) - 7l(k)74k)) (4.9) 



where z, a unit vector orthogonal to the x, y plane. As in I, J/ measures the current 
carried by motion of the quasiparticles and quasiholes — notice the A prefactor, indicating 
that these excitations have opposite charges. The operator J/J creates a quasiparticle- 
quasihole pair: throughout this section we will work in the approximation in which J// = 0. 
We will also neglect all mixing between the <^7±7±^ and <^7±7q:^ distribution functions. 
These approximations simply amount to neglecting all but the a^^ component of a — the 



other components were already computed satisfactorily in perturbation theory in Section |T1 



We can now use the standard equation of motion analysis ||2^ to write down the colli- 
sionless transport equations for the excitations. As a first step, we define the distribution 
functions 

/,(k,t) = (7l(k,t)7A(k,t)). (4.10) 

In equilibrium in the absence of external perturbations, these are related to the Fermi 
function defined in (|3l 
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U{k,t) = f{k) 

f_{k,t) = f{-k) = l-f{k) (4.11) 

Then to first order in a, in the presence of an external electric field E, we find the simple 
equations 

At this order, in the approximations discussed above, the only effect of the interactions has 
been to replace E by E^. It is a simple matter to solve (|4.12|) in linear response. First we 
parametrize the change in fx from its equilibrium value by 

/A(k, uj) = 27T6{uj)f{Xk) + Xqek ■ ET{uj)g{k, u), (4.13) 

where we have performed a Fourier transform in time to frequencies, u, and introduced 
the unknown function g{k,uj). Notice that the change in the distribution functions has an 
opposite sign for quasiparticles and quasiholes. Thus the quasiparticles and quasiholes move 
in opposite directions: as they have opposite charges, their electrical currents are equal, 
while their net momenta have opposite signs. Inserting (|4.13| ) into (|4.12| ), we can obtain a 
simple solution for the function g 

Inserting this result into ( [4.10| ) and (|4.8|), we obtain the conductivity 

which is precisely our earlier result (|3.13|) , in the approximation M{T) = and J^{a) = ln(2). 
In this approach, it is now evident that the leading factor of 2 comes from the sum over A: 
i.e. the quasiparticles and quasiholes contribute equally to the total current. 



B. Collision-dominated transport 



We now include collision terms on the right hand side of (|4.12|) . As in I, while the terms 
can be formally derived using the methods of Ref , it is easier to determine by a simple 
application of Fermi's golden rule. In the latter manner we obtained 



'| + ,eE..-l/.(k,«) 



dk) 

(Pki (fq 



{2n)6{k + k, - |k + q| - |ki -q|)< 



(27r)2 (27r)2 

\T,\' {/,(k, t)/-A(-ki + q, t) [1 - Mk + q, t)] [1 - f-x{-kut)] 
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- [1 - /A(k,t)][l - /_,(-ki +q,t)]/,(k + q,t)/_,(-ki,t)} 
+ ^ {/,(k, t)/,(ki, t) [1 - /,(k + q, t)] [1 - /,(ki - q, t)] 

- [1 - /A(k,t)][l - /A(ki,t)]/,(k + q,t)/,(ki - q,t)}|. (4.16) 



where 



Ti = T+__+(k, -ki + q, q) - T_+_+(-ki + q, k, k + ki) 

- r+_+_(k, -ki + q, -k - ki) + r_++_(-ki + q, k, -q) 
^2(k, ki, q) = T++++(k, ki, q) - T++++(ki, k, k + q - ki) 

- T++++(k, ki, ki - k - q) + T++++(ki, k, -q) (4.17) 

The terms proportional to |Tip represent collisions between oppositely charged particles, 
while those proportional to are collisions between like charges. There are also processes 
a particle-hole pair is created: as in I, these can be dropped because they have vanishing 
phase space upon imposition of the energy conservation constraint with dispersion Sk = k 

We now proceed to the linearization of ( [4.16| ) by inserting the parametrization ( [4.1 3|) 
and find 

-iujkg{k,uj) H = 



(27r)2(2vr)2^^^)^^^ + ^i~l^ + 'll~l^^"^l) 

g(|ki-q| + |k+q|)/T ^ g/ci/T)(|J.^|2 ^ 1^2172) 

eki/T + i)(e|ki-q|/T + i)(elk+q|/T + 1) 

(|ki-q| + |k+q|)/T ^ e'=/^)(|Ti|2 - IT2I72) 
gfc/T + l)(e|kl-q|/T + l)(e|k+q|/T + 1) 
elk+q|/T ^ g(fc+fci)/T)(|y^|2 _ 1^2172) \ 

(efc/T + i)(efci/r+i)(elk+q|/r + i) j 

elki-q|/T + g(fe+fei)/T)(|j.^|2 ^ |T2|72)~ 

(efc/T + i)(e*=i/T + l)(elki-q|/T + i) 



+ 





(4.18) 



The final ingredient necessary before assembling all the results of this paper is now a solution 
of ( 4.18| ) for the distribution function g{k, uj). An analytic solution is clearly out of the ques- 



tion and, we obtained the solution numerically. However, before attempting any numerical 
solution, the singular nature of the scattering cross-sections |Tip and IT2P makes it essential 
to obtain an analytic understanding of the nature of any potential infrared divergences aris- 
ing from the momentum integration. We will show below that all the potentially infrared 
divergent terms in the collision integral do indeed vanish. However just this knowledge is 
not sufficient to allow one to proceed with the numerical integration. One does not want 
numerically integrate an integrand with potentially divergent terms which will eventually 
average out to a finite answer, as such a procedure is clearly prone to large systematic nu- 
merical errors; rather one wants an integrand which is explicitly free of singularities. On the 
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other hand, the integrand is dependent upon the unknown function g, so exphcitly removing 
the singularities is difficult. These issues are clearly crucial, but also vexing, and stymied 
the author for some time. Eventually, a remarkably simple and useful parametrization of 
the collision term was found. As the analytic discussion of the cancellation of singularities 
is also simplest in this new parametrization, we will ffist introduce it, and then discuss the 
cancellation. 

The key is to carry out the q integration in the collision term of ( [4.18[ ) ffist, and to notice 
that the delta function associated with energy conservation 

5(fc + A;i-|k + q|-|ki-q|) (4.19) 

defines an ellipse in the q plane with foci at — k and ki. It therefore pays to parametrize q 
in an elliptic co-ordinate system p5|. We introduce elliptic co-ordinates 

0<fi<oo , -n <e <7r (4.20) 

and express q in terms of fi and 6 by 

Q = ^'-^ + ^JSl±I^ (eM+^^ + e---^) . (4.21) 

Here k and ki are regarded as fixed external parameters, and ki will be subsequently 
integrated over. Then it can be easily verified that the ellipse represented by (|4.19| ) is 
given simply by 

cosh/i = - 4.22 

|k + ki| 

The integral over q in ( [4. 181 ) is a two-dimensional integral over /i and 6, and the integral 
over n can now be trivially carried out. This is summarized by the general formula 



72 

I -^{2n)5ik + A;i - |k + q| - |ki - q|)F(q, k, k^ 



r2Tr dQ roo Ik + kx| / \ 

/ — / dfx (cosh^/i — cos^ 6') (5(/i; + /ci — |k + ki| cosh/i)F(Q, K, i^Ti) 

JO 27r JO 4 ^ ^ 

27r Ik + kil fcosh^ u - cos^ O) 

>-F(Q,K.K,), (4.23) 

where F is an arbitrary function, and in the last equation we use the value of /i given in 
(g) and Q given in (g) 



We are now prepared to discuss the handling of the infrared divergences in ( [4.18|) . There 
are two distinct types of potential divergences and we will consider them in turn. 

• The ffist singularity is associated with the divergence in the forward scattering cross- 
section. By an examination of ([4.17|) we see that there is a divergence in |Tip 
when the momentum transfer q approaches zero, and also a l/(q — ki + k)^ divergence 
in \T2\^ when the momentum transfer q approaches q = ki — k; the latter divergence 
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corresponds to forward scattering with an exchange of particles. (Note that there is 
no divergence when the momentum transfer k + ki vanishes, as is easily seen by an 
examination of the A labels in ( [4.17| ), combined with the expression ( |4.6| )). However, 



we then notice from the right hand side of (|4.18| ) that the coefficient of |Tip in the 



collision term vanishes at q = 0; the term in the co-efficient linear in q will vanish after 
an appropriate angular integral over q, and the remaining combination is explicitly 
finite at q = 0. Similar considerations apply to the the coefficient of IT2P, which 
vanishes at q = ki — k. The physics of this cancellation is simple and quite familiar. 
Its origin is the "(1 — cos (scattering angle))" phase space factor that appears in the 
expression for transport relaxation rates this vanishes at zero scattering angle, 
and cancels out the divergence in the scattering cross-section. 

Now there remains the issue of how to numerically treat this potential divergence. It is 
clear that the angular integral around q = and q = ki — k should be taken with great 
care, or we will be left with a spurious divergence. Here our elliptical parametrization 
( [4.21|) turns out to be helpful. Notice that, as must be the case, the points q = and 



q = ki — k satisfy the delta function in ( [4.19|) , and therefore lie on the ellipse defined 



by ( [4.221) . These points lie at = and 6 = 6q -\- n where do is determined by the 



equation 



This is a quadratic equation for e*^", and therefore has two solutions which are easily 
determined. For fi give by ( [4.22| ), one of the solution will have unit modulus and 
therefore determines a real 6q: this is guaranteed by the fact that q = lies on the 
ellipse defined by ( [4.221 ). The angular averaging is carried out by the 6 integral in 
( |4.23|) , and we now simply have to insure that the vicinities of the points 6 = 60 (and 
= 60 + Tc) are treated symmetrically. More precisely, we always pair the integrands 
at + ^ and 6q — 6 { and also at ^0 + ^i" + ^ and 6q + n — 9 ) so that the combined 
integrand is explicitly finite at the potentially singular points. 

The second potential singularity arises purely from phase space considerations, and 
is a curious property of energy-conserving scattering of particles in two dimensions 
with the linear dispersion Sk = k. We are considering the scattering of two incoming 
particles with momenta k and ki. Consider the case where these two momenta are 
collinear and point in the same direction i.e. when k ■ ki = kki. We are eventually 
integrating over ki in the collision term of ( 4.1^ ), and this condition defines a semi- 



infinite line in the two-dimensional ki space. Along this line, we see from ( |4.22| ) that 
/i = 0. Now the phase space factor 1/ sinh/x in ( |4.23|) looks dangerous (generically the 



q integral in (|4.23[ ) will be non-singular, and will not effect the considerations here). 
This infinite factor is present along the entire semi-infinite line in the ki plane: if 
we pick a point a distance y from this line (almost collinear particles) we find that 
the 1/ sinh/i factor diverges as 1/y. This implies a logarithmic divergence in the ki 
integral. Fortunately, the specific form of Ti, T2 under consideration here comes to the 
rescue. It is easy to show from (|4.6| , [4.17D that these cross-sections vanish identically 
for the case of the scattering of particles moving in the same direction. 
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The remaining steps in the solution of the integral equation were relatively straight- 
forward. It is easy to show that g{k^uj) ~ e~^l'^ for large /c, and so it is an excellent 
approximation to simply truncate the infinite range of k integration at some large positive 
k. Over this finite range, we expanded gik^uS)^!'^ over the space of Chebyshev polynomials. 
We determined the action of the kernel of the integral equation on each Chebyshev poly- 
nomials, and decomposed the resultant also over Chebyshev polynomials. This defined a 
discrete matrix, which is the transcription of the kernel to the discrete basis space of Cheby- 
shev polynomials. Solving the integral equation then simply became a matter of inverting 
the matrix. Notice that as the kernel does not explicitly depend upon frequency, the same 
matrix could be used at any frequency. 

We now present the results of the above numerical analysis. A simple dimensional anal- 
ysis of ( [4.18| ) shows that the entire T and a dependence of the integral equation can be 
scaled out, and the resulting scaled equation is parameter free. In particular the solution of 
( [4.18| ) has the form 



g{k,uj) 



1 



:G 



ijj 



a- 



T3 XT' a^T 



(4.25) 



where G is a universal function we have determined numerically. Inserting the form ( [4.25| ) 
into ( |4.13| ) and ( [4.8|) we see that the conductivity ( |3. 131 , ^1151) is now replaced by 



- ilk) 



(4.26) 



with the universal scaling function S^p given by 



k^dkGik,u) 



(4.27) 



This scaling function describes the broadening of the delta function in the real part of ( [4.151) : 
the large frequency behavior of ( [4.18| ) shows that it has the same total spectral weight as 
(H): 



-Re 



TT 



ln(2) 



(4.28) 



We show the numerical solutions for the universal functions G and S^p in Figs |^ and 
respectively. The value of SJ^p(O) determines the d.c. conductivity 



a^:,{u/T 0) 



gV \ 0.437 



(4.29) 



This result and ( |3.17| ) and ( |3.28| ) completely specify the limiting forms of the conductivity 
a. 
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V. SYNTHESIS 



We will now finally turn to a discussion of the properties of the physical conductivity 
cr. The several disparate ingredients necessary for its computation have all been assembled, 
and let us now summarize them. We began with a model with two species of anyons il'i,2- 
We defined their irreducible conductivities ai^2 and argued that the physical conductivity a 
was related to them by ( p.l3|) . The iIj2 anyons were non-critical and the value of 5^2 is given 
simply by ( |2.14| ). Sections Q and |IV| were devoted to studying the value of the critical ai. 
Section |ITI| considered the perturbative contribution which in the limit u/T oo reduced 
to the coherent transport of externally created particle-hole pairs. The Hall conductivity 
was given entirely by such a contribution: a^^j + o-^J + d'j^y where the results are given in 
( 3.191 ), (|3.23| ) and ( p. 24 ) respectively. The longitudinal conductivity consists of the coherent 



contribution cr™^ given in ( |3.15| ) , and an incoherent quasiparticle contribution ct^p which was 
computed in Section the latter satisfies the scaling form ([4.26|) and the scaling function 



E^P was plotted in Fig |. 

The above discussion can be restated as follows. Let us define the dimensionless functions 
A 5 by 

B{u) = -1 + 4-2 KH^) + + Sfgn^)) (5.1) 

Then, the physical conductivities are given by 



h / a^A^iuj) + (l + a5(cj))2 



^^^^""^ - ) a^A^{uj) + {l + aBiuj)y ^^'^^ 

We emphasize that A and B are complex, universal functions of uj/T, and it is necessary to 
include both their real and imaginary parts in the above expressions. 

In evaluating our final results for a^x, an important complicating factor becomes appar- 
ent. While we are ultimately interested in a of order unity, the terms entering (|5.2| ) are of 
many differing orders in a for small a. As a result, certain processes become unphysically 
dominant at small a, while in the real system we expect a more even-handed competition. 
Further the terms also vary on differing scales of u. In particular, the quasiparticle/quasihole 
conductivity a^^ is of overall order 1/a^ and varies on a frequency scale uj ~ a^T; for small 
uj it is the most dominant term as a — 0. Of the other terms, is order unity, while 
a^J^'^'^^ are of order a, and all vary on a scale uj T] these become important for large uj 
as a ^ 0. Because of these complications, we believe more significance should be attached 
to the form of our results for the individual conductivities in Section |TI| and The final 
results, when these distinct components are combined, are expected to be less reliable for 
moderate a. 

We begin by describing the formal structure of the small a expansion for a. It is necessary 
to consider two distinct frequency regimes separately: 
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(A) uj of order a^T, but with uj <^T: In this ct^p is the dominant term, and examination of 
the small a limit of (|5.2|) gives 



T.Tx{uj/a^T) 



+ ... 



(5.3) 



We show a plot of the result above for Re [(Ta;2,(ci;)] in Fig ^. Notice that this not simply 
the inverse of the real plot in Fig ^, as the first equation in ( ^.31) involves the inverse of the 
comj^/ex function S^p. The function Re l/E^P(a;) approaches a finite number in both the 
limits tl; — > (where it is simply the inverse of ( [4.29|) and tl; oo: 

1 



Re 
Re 



l^xx \UJ 
1 



0) 



^xx \ 



2.29 
3.47 



'~{iU oo) 

In a similar manner, we can also get the Hall conductivity in this regime: 

1 



(5.4) 



crxy[^) 



a 



+ 



(5.5) 



Surprisingly, this result for a^y has the opposite sign from the Hall conductivity of the 
quantized Hall phase connected to the critical point under consideration here: we strongly 
suspect this feature is an artifact of the absence scattering from interactions other than the 
Chern Simons term: a model in which a^P is not so singular in the small a limit should have 
a Hall conductivity of the "correct" sign. 

(B) UJ of order or larger than T: Now the result is constructed simply by inserting the per- 
turbative results of Section |Tl| for a in (p.2|). The formal small a expansion is 



(TxxK,^^) 
C^xyi.^) 



2h 



+ CTxy + 



q^e^ 



+ 



(5.6) 



and the ingredients determining the coefficients were plotted in Figs |T] and ^ 

An alternative approach to gaining some understanding of the predictions of (|5.2| ) is to 
simply evaluate the full expression at a given value of a. This has been done in Figs ^ and |^, 
which show the results of directly evaluating (p.2|) at face value at a = 0.3. 



VI. CONCLUSIONS 



The purpose of this paper was to describe the dynamic transport properties of the the 
simplest possible quantum critical point in a quantum Hall system. In particular, motivated 
by experimental results discussed in Section |, we chose a model quantum phase transition 
which had non-zero interactions in the critical fixed point theory. The critical properties of 
such theories have a universal dependence of all observables on uj/T, and, as pointed in I, 



22 



the dependence of the conductivity on u /T at two dimensional quantum critical points is 
especially important. The model we chose, S in (|2.1|) , is equivalent to a model of free anyons 
and anti-anyons with a relativistic Dirac spectrum. This model is properly considered as 
an interacting theory for the following reason: two anyonic/anti-anyonic excitatios have a 
non-trivial T matrix for scattering in which the momenta of the ingoing and outgoing 
particles can be arbitrary, and are subject only to the constraints of conservation of total 
energy and momentum (it is clear that this fails only for purely bosonic or fermionic statis- 
tics, when the ingoing and outgoing momenta are the same). Such scattering is the main 
ingredient in establishing local equilibrium, and is essential to the description of transport 
in the regime -C T. 

Our results for the non-zero temperature transport properties were summarized in Sec- 
tion [V]. It is evident from the discussion there that, while the general scaling forms for 
the conductivities are very similar to those discussed in I for the boson superfiuid-insulator 
transition, the explicit results for the scaling functions themselves are quite different. For 
instance, we found in I that the longitudinal conductivity a^xi.^) was a decreasing function 
of increasing uj near = 0. In the present quantum Hall problem, we found that for small 
a, the longitudinal conductivity satisfied ( |5.3|) , and is therefore seen to be an increasing 
function of uj near u; = 0, as shown in Fig ^. This property was subsequently responsible 
for a peak in a^x near uj T shown in Fig ^ 

As they stand, our model calculations can probably not be directly applied to microwave 
conductivity measurements in the quantum Hall system [|^. The most important defect 
that must be repaired is the absence of Coulomb interactions among the charged excitations. 
A simple power counting argument shows that the Coulomb interaction is a marginal per- 
turbation at the fixed point considered here. A theory for the consequences of this marginal 
perturbation is not available at present, and appears to us to be an important subject for 
future research. It is interesting to note that since the Coulomb interaction is marginal at 
tree level, it must be characterized by a dimensionless coupling constant. It is easy to see 
that this coupling constant is 

#=f#l--^ (6-1) 
nv \nc J V 137 

where v is the velocity of the Dirac fermions, c is the velocity of light, and the quantity in 
the parenthesis is the dimensionless fine structure constant. For v, it seems reasonable to 
use a typical drift velocity of electrons in the two-dimensional electron gas ~ 10^ — lO^m/s, 
and we see then that this dimensionless coupling can be quite large. The Dirac fermion 
model (at a = 0) also describes quasiparticles in a d-wave superconductor, and so the 
dimensionless coupling characterizing the Coulomb interactions can also be quite large for 
the high temperature superconductors. 
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APPENDIX A: LARGE N EXPANSION 



Our study of transport near the superfluid-insulator quantum critical point in I, and of 
the quantum Hall critical point in the present paper, has been carried out using transport 
equations in which a four-point coupling between the charge carriers is used as an expansion 
parameter. In this appendix we will outline an alternative approach to the same problems, 
the large expansion. We shall not provide an explicit numerical solution of the large N 
transport equations here, but it will evident that the structure of the results is very similar 
to those already obtained. 

We will begin with the simpler case of the superfluid-insulator transition in a model 
of 0(2) quantum rotors which was discussed in I. We generalize this to a model of 0{N) 
quantum rotors, and the quantum field theory describing the critical point has the action 

Sr = fj^ dr j£x[^ [idM' + c'iVAaf + ml<Pl] + . (Al) 

Here (pa is a A^-component field and the action has 0{N) symmetry (the 0{N) index a is 
implicitly summer over). The spatial and temporal gradient terms are both second order, 
so the action has a "Lorentz" invariance with c the velocity of light, and as a result the 
dynamic critical exponent z = 1; we will henceforth use units in which c = 1. The quartic 
non-linearity, uq, has been scaled by a factor of to allow the interesting large N limit. 
The bare "mass" term, is cutoff dependent, and must be adjusted order by order in uq so 
that the system is at its quantum critical point at T = 0. The superfluid-insulator transition 
of boson Hubbard models is described by the case N = 2, where \1/ = 0i + z</)2 is the usual 
complex superfiuid order parameter. The case N = 3 which applies to quantum-critical 
points in quantum antiferromagnets p2|j27[ . 



First we review the large expansion for the order parameter susceptibility, which is 
the propagator, G, of the 0a field. This was extensively discussed in Ref At momentum 
k, and Matsubara frequency Un, the N = oo result for this propagator is 

G{Ktu^-) = ,.^J^Q.j. , (A2) 
where is a pure (universal) number: 

e = 2ln(^^] (A3) 



The corrections are obtained by decoupling the quartic term in ( [AID by a Hubbard- 



Stratanovish field A (not to be confused with the charge subscript A used in Section |V|, and 
also below). There is then a three-point coupling ~ 0^A in the action. The propagator of 
the A field is given by the inverse of the bubble graph of the (pa field. In the low-energy, 
long-distance limit relevant for the critical properties, this propagator becomes independent 
of the microscopic value of uq, and equals 

2 1 

iVH(q,^eJ ^ ^> 
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where 

f <Pk 

n(q, ien) J j^^G(k, iuJn)G{k + q, iuj^ + ien) (A5) 

Notice that there is a factor of associated with each A propagator: this comes from the 
sum over the N components of the (pa field in the bubble graph. So the 1/N expansion for 
the susceptibility is generated by expanding in the number of A propagators. 

Now we turn to transport properties. The value of the conductivity at = oo was 
discussed in Section II of I. The structure was very similar to that obtained in an expansion 
in e = 3 — d: there was a quasiparticle/quasihole delta function at = 0, and a continuum 
contribution at frequencies greater than u = 2BT. The main purpose of this appendix is to 
determine the manner in which the u = delta function broadens out in the expansion. 
Naturally, this will be done by formulating the transport equations of the quasiparticles and 
quasiholes in the large N limit. 

The general transport equations for Sr were discussed in Appendix B of I in an 
expansion in uq. Following that treatment, we consider here the linear response to an 
"electric" field, E, which couples only to the 0i 2 components i.e. the corresponding electric 
potential couples to a generator of 0{N) which rotates the system in the 1,2 plane. Then 
it was argued that the distribution function of excitations in the remaining components 
remained unchanged to linear order in the field — all corrections are at least quadratic in E. 
In other words, if, as in I, we make the mode expansion in terms of particle creation and 
annihilation operators 

0a(x, ^) = / (0 {M^^ ty'"'' + ^I(k, t)e-'''-') , (A6) 
with the dispersion 

= + Q^T^f'\ (A7) 

then, we have 

/,(k,t) = (6l(k,t)6,(k,t)) = n{eu) = for a > 2. (A8) 

This has very important consequences for the large A^ expansion. Recall that the inverse 
propagator of the A field was given by the sum over the bubble graphs of the A^ components 
of the (pa- In the large A^ limit, this sum is dominated by the a > 2 components, and so also 
remains unaffected to linear order in E. So we can assume that the propagator for A is still 
given by (|A4|) . 

Now the quantum transport equations can be easily written down by noticing that the 



present problem has a structure closely analogous to the electron-phonon problem p6| , |28 
the A field plays a role similar to the phonons, and the 0^A coupling is analogous to the 
electron-phonon coupling. The fact that the A propagator is not modified by E is equivalent 
to the statement that "phonon drag" effects are negligible in the large A^ limit. Indeed, the 
particular form of the propagator in ( [A4D translates into a special phonon density of states, 
and the transport equations can then be read directly off Ref p8|. Let us define 
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V2 



(A9) 



Then fx(k,t) = (b\{'k,t)bx{'k,t)'^, with A = ± satisfies the transport equation 
f d , ^-^ 5 \ , „ , 2 dVt f (fq ( 1 \ 



X ■ 



{2'K)5{ek - ^ik+qi - fi) 



k+q| 



/A(k,t)(l + /,(k + q,t))(l + n(fi)) 



+ 



-/A(k + q,t)(l + /A(k,t)H(]) 
(27r)5(efc - 5|k+q| + Vt) 



4£fc£|k+q| 



^(k,t)(l + ^(k + q,^))r^(^]) 



+ 



-/A(k + q,t)(l + /,(k,t))(l + n(fi)) 
(27r)5(£fc + £|_k+q| - 



4£fc£^|-k+q| 



/A(k,t)/_A(-k + q,t)(l + n(fi)) 



- (1 + /_A(-k + q, + /A(k, 



(AlO) 



where Q is the charge of the ha quanta. Now this equation can be hnearized hke ( ^.16 ), and 



then solved using the numerical methods discussed in Section IV. We will not present the 



numerical solution here, but an essential fact is already clear from the structure of (|A10[ ). 
The u; = delta function in the conductivity, found at = oo, gets broadened out to a 
width uo ~ T /N , and the peak conductivity at = is of order A^. 

The above large A^ analysis can obviously be carried out at any 1 < < 3. It is useful to 
do this and then take the limit of small 3 — c?. We find that the equation so obtained agrees 
precisely with the large A^ limit of the small ?> — d equation discussed in Appendix B of I 
for general A^. Thus the N ^ oo and 3 — ^ limits of the quantum transport equations 
commute, and this is an important check of the consistency of our analysis. 

Our results can be easily extended to the large A^ limit of either bosonic |18| or 
fermionic |]Ty| theories with Chern-Simons terms. For instance, an alternative formulation 
of the theory of anyons considered in this paper is to represent them as relativistic bosons 
with attached statistical fiux tubes: this is represented by a generalization of the action Sr 
as follows [0. Consider the case with A^ = 2M even, parametrize the 2A^ real fields (pa in 
terms of M complex fields ^ 



+ i(j)2a)/V2 a 



M 



(All) 



Now couple a U{1) Chern Simons field to the diagonal U{1) charge 

Sra = dT J d^x[\{dr - iar)^a\^ + C^\{Vx - ia^)^a\^ + mlJ^al 



(A12) 
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Notice that the effective statistcs parameter of the anyons is a/M: this is necessarily very 
close of the bosonic limit for large M. 

The large M analysis of S^a is very similar to that of S: in addition to the field A, we 
now have a gauge field whose propagator is also given by a boson current polarization 
bubble, and is of order 1/M |18[. By arguments very similar to those for A, we can conclude 
that the deviations in the propagator due to the external electric field can be neglected 
i.e. "photon drag" is negligible in the large M limit. The resulting quantum transport 
equation is then very similar ( [AlOj ): along with the terms already displayed, we have an 
additional parallel set with Im(l/n) replaced by the spectral density of the propagator 
times current matrix elements of the \l/ field. 

Finally, similar considerations apply to the large M generalization of the fermionic Chern- 
Simons theory, S, of this paper. We replicate the fermionic terms in iS M times, rescale 
a —>■ a/M, and then take the large M limit as discussed above for the bosonic theories. 



APPENDIX B: COMPUTATIONS FOR THE HALL CONDUCTIVITY 



This appendix shows how the result (|3.20| ) is obtained from ( p. 191 ) . First it is clear after 

de 



using (pl\ ) and rescaling the integrand by T that 



^Si^/T - 0) = ^ ln(2) / ^ tanh(e/2) 



h 



Let us define 



F{y)= ^tanh(x/2). 



. . tanh(y/2) ln(?//2) /""^ , ln(x) sinh(x) 

F{y) = ,9. + / dx- 



y 2cosh^(y/2) Jy/2 cosh^(a;) 

Now we can explicitly take the y limit to get 

^ , 1 , , s ln(2) + 1 , ln(x) sinh(x) 

F{y-.0) = -\n{l/y) + ^ +/ dx ^ > > 

I Z Jo cosh [x] 



(Bl) 



(B2) 



and it is clear that to understand the small a limit of ([BID, need the form of F{y — > 0). 
By integrating twice by parts, it is easy to show that 



(B3) 



(B4) 



Then, we evaluate the last integral by first expanding the integrand as an infinite series in 



-2x. 



°° ^ ln(x) sinh(x) 



cosh (x) 



-4 / dx\n{x) 
Jo 



E(-i: 



.n=l 



2^(-l)"n(7 + ln(2n)) 



n=l 

oo 



2 ^(-1)" ((7 + ln(2))ni+^ ' ^ 



n=l 



ds 



-n 
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6(7 + ln(2))C(-l 
71n(2) 



7 
2 



6 



ds 

-6C'(-1). 



2- ((2^ 



-i)c(-i-^; 



(B5) 



Notice that the intermediate steps involve badly divergent series in n, and so the validity 
of the final result is not a priori obvious. To check it, we numerically evaluated the left 
hand side and verified that it was equal to the final right hand side to a very high accuracy. 
Combining (0), (|B2D, (0) and we can now obtain ( CT) - 
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FIG. 1. Real (full line) and imaginary (dashed line) parts of the perturbative result ( |3.15| ) for 
the conductivity a^'^{uj) as a function of lu/T evaluated at a = 0.3. The conductivity is measured 
in units of q^e^ jh. The singularity is at w = 2M{T), and is artifact of the absence of damping 
at this order. Note, as in ( |3.12 ), the conductivity axx contains another quasiparticle/quasihole 
contribution which is a delta function in the perturbative approach. 
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FIG. 2. Real and imaginary parts of the perturbative results for the Hall conductivity ( ^.191) , 



( 3.23| ) and ( 3.24| ) for axy{uj) = a^y + ^xy as a function of the uj/T evaluated at a = 0.3. The 

conductivity is measured in units of q^e^ /h. Again the spurious singularities are at a; = 2M{T). 
Note that unlike, (Txx, the Hall conductivity is adequately described by the perturbation theory. 
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FIG. 3. Real and imaginary parts of the solution for the universal function /c^G/T^, introduced 
in ( [4.25| ), as a function of k/T for a few values oiu/a^T. The combination k'^G/T'^ is that appearing 
in the integrand of the integral ( [4.27 ) for the conductivity. 
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FIG. 6. Real and imaginary parts of the axx as a function of co/a'^T, plotted at a = 0.3. The 
conductivity is measured in units of q^e^ /h. The sing ularities at w = 2M{T) {uj/a^T ^ 9.2) will 
be rounded out at higher orders. 
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FIG. 7. Real and imaginary parts of the function of w/a^T, plotted at a = 0.3. Again, 

the conductivity is measured in units of c^e^ jh and the sing ularities at w = 2M(r) {uj/a^T ^ 9.2) 
will be rounded out at higher orders. 
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